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We report measurements of B to pseudoscalar-pseudoscalar decays with at least one ij meson in 
the final state using 140 fb _1 of data collected by the Belle detector at KEKB e + e~ collider. We 
observe the decay B + — > rjn + and find evidence of B + —> nK + ; the measured branching fractions are 
B(B + -> r/ir + ) = (4.8io.?(stat)±0.3(sys)) x 1(T 6 and B(B + -> r)K+) = (2.1 ±0.6(stat) ±0.2(sys)) x 
10 -6 . Their corresponding CP violating asymmetries are measured to be 0.07±0.15(stat)±0.03(sys) 
for ■i]tt ± and —0.49 ± 0.31 (stat) ± 0.07(sys) for r/K ± . No significant signals are found for neutral 
B — > r\h decays. We report the following upper limits on branching fractions at the 90% confidence 
level: J3(B° -> V K°) < 2.0 x 10~ 6 , J3(B° -> r?7r°) < 2.5 x 10~ 6 and B(B° ->■ rjrj) < 2.0 x 10" 6 . 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Charmless B decays provide a rich sample to understand B decay dynamics and to search for CP violation. An 
unexpectedly large B — > rjK branching fraction 0,0 has stimulated much theoretical interest. It was suggested even 
before the rjK measurement that two b — ► s penguin amplitudes are constructive in B — » rj' K decays but destructive 
in B — > rjK The situation is reversed for B — > rjK* and B — ► rjK* decays. Experimental results have more 
or less confirmed this picture; however, precise measurements of branching fractions are needed to quantitatively 
understand the contribution of each diagram. It was also pointed out that in the rjK mode the suppressed penguin 
amplitudes may interfere with the CKM suppressed b — ► u (tree) amplitude and result in direct CP violation Q. 
The penguin-tree interference may also be large in B + — > i]'tt + 5] and B + — > 7]ir + decays; however, theoretical 
expectations for the partial rate asymmetry (Aqp) can be either positive or negative Recently, the BaBar 

Collaboration has reported large negative Aqp values in both rjK + and i]ir + , which are ~ 2a away from zero 0. 
However, more data are needed to verify these large CP violating asymmetries. Furthermore, branching fractions 
and partial rate asymmetries in charmless B decays can be used to understand the tree and penguin contributions 
and provide constraints on the third unitarity triangle angle <f>3 ||. 

In this paper, we report measurements of branching fractions and partial rate asymmetries for B — » rjh decays, 
where h could be a K, n or 77 meson. The partial rate asymmetry is measured for the charged B decays and defined 
to be: 

_N(B~ ->rih-)-N(B + rjh + ) 
CP ~ N(B~ -> r]h-) + N{B+ -> rjh+) ' 

where N(B~) is the yield for the B~ — > rjh~ decay and N(B + ) denotes that of the charge conjugate mode. The data 
sample consists of 152 million BB pairs (140 fb _1 ) collected with the Belle detector at the KEKB e + e~ asymmetric- 
energy (3.5 on 8 GeV) collider || operating at the T(45) resonance. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three-layer silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cercnkov counters (ACC), a barrel- like 
arrangement of timc-of-flight scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of 
CsI(Tl) crystals located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons (KLM). The detector 
is described in detail elsewhere [inj . 

Two ij decay channels are considered in this analysis: rj — > 77 (?y 77 ) and rj — > n + ir~ir° (r/^). In the ?7 77 reconstruc- 
tion, each photon is required to have a minimum laboratory energy of 50 MeV and the energy asymmetry, defined as 
the absolute value of the energy difference in the laboratory frame between the two photons divided by their energy 
sum, must be less than 0.9. Furthermore, we remove 77 candidates if either one of the daughter photons can pair with 
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any other photon to form a ir° candidate. Candidate 773^ mesons are reconstructed by combining a 7r° with a pair of 
oppositely charged tracks that originate from the interaction point (IP). We make the following requirements for the 
invariant mass on the ij candidates: 516 MeV/c 2 < M 77 < 569 MeV/c 2 for ?7 77 and 539 MeV/c 2 < M 37r < 556 MeV/c 2 
for 77 37r . After the selection of each candidate, an rj mass constraint is implemented by readjusting the momenta of 
the daughter particles. 

Candidate neutral pions are selected by requiring the two-photon invariant mass to be in the mass window between 
115 MeV/c 2 and 152 MeV/c 2 . The momentum of each photon is then readjusted to constrain the mass of the 
photon pair to be the nominal ir° mass. To reduce the low energy photon background, each photon is required 
to have a minimum energy of 50 MeV and the 7r° momentum must be above 250 MeV/c in the laboratory frame. 
Charged tracks are required to come from the IP. Charged kaons and pions that form B candidates with 77 mesons 
are identified by combining information from the CDC (dE/dx), the TOF and the ACC to form a K(n) likelihood 
Lk{L-k)- Discrimination between kaons and pions is achieved through the likelihood ratio Lk/{L-x + Lk)- Charged 
tracks with likelihood ratios greater than 0.6 are regarded as kaons, and less than 0.4 as pions. The typical kaon 
and pion identification efficiencies for 2.5 GeV/c tracks are (85.0 ± 0.2)% and (89.3 ± 0.2)%, respectively. With the 
same track momentum, the rate for pions to be misidentified as kaons is (7.3 ± 0.2)% while the rate for kaons to 
be misidentified as pions is (10.6 ± 0.2)%. Furthermore, charged tracks that are positively identified as electrons or 
muons are rejected. Kg candidates are reconstructed from pairs of oppositely charged tracks with invariant mass 
(M^) between 480 to 516 MeV/c 2 . Each candidate must have a displaced vertex with a flight direction consistent 
with a Kg originating from the IP. 

Candidate B mesons are identified using the beam constrained mass, Afbc = \J E 2 c&m — P B , and the energy differ- 
ence, AE = Eb — Ebcam, where -Ebeam is the run-dependent beam energy in the T(45) rest frame and is determined 
from B — ► D^n events, and Pb and Eb are the momentum and energy of the B candidate in the T(45 l ) rest frame. 
The resolutions on Mbc and AE are around 3 MeV/c 2 and ~ 20-30 MeV, respectively. Events with Mb c > 5.2 
GeV/c 2 and \ AE\ < 0.3 GeV are selected for the analysis. 

The dominant background comes from the e + e~ — > qq continuum, where q = u, d, s or c. To distinguish signal from 
the jet-like continuum background, event shape variables and the B flavor tagging information are employed. We form 
a Fisher discriminant IllJ from seven variables that quantify event topology. The Fisher variables include the angle 0t 
between the thrust axis |l2j of the B candidate and the thrust axis of the rest of the event, five modified Fox- Wolfram 
moments 01, and a measure of the momentum transverse to the event thrust axis (S±) [l4|- The probability density 
functions (PDF) for this discriminant and cos 9b, where 8b is the angle between the B flight direction and the beam 
direction in the T(45) rest frame, are obtained using events in signal Monte Carlo (MC) and data with M^ c < 5.26 
GeV/c 2 for signal and qq background, respectively. These two variables are then combined to form a likelihood ratio 
1Z = £ s /(£ s + C-qq), where £ s (qq) is the product of signal (qq) probability densities. 

Additional background discrimination is provided by the quality of thei? flavor tagging of the accompanying B 
meson. We use the standard Belle B tagging package |l5j, which gives two outputs: a discrete variable (q) indicating 
the B flavor and a dilution factor (r) ranging from zero for no flavor information to unity for unambiguous flavor 
assignment. We divide the data into six r regions. Continuum suppression is achieved by applying a mode dependent 

requirement on 1Z for events in each r region based on N^ p / Ng Xp + N°* p , where A r ° xp is the expected signal from 

MC and N^ p denotes the number of background events estimated from data. This 1Z requirement retains 58-86% 
of the signal while reducing 96-82% of the background. From MC all other backgrounds are found to be negligible 
except for the rjK + <-> t]tt + reflection, due to K + «-> 7r + misidentification, and the r\K *(892)(r;p(770)) feed-down to 
the r]K(r]ir) modes. We include these two components in the fit used to extract the signal. 

The signal yields and branching fractions are obtained using an extended unbinned maximum-likelihood (ML) fit 
with input variables Mb c and AE. The likelihood is defined as: 

N 

C = exp (- J2 N j) II [ E N j P j( M ^> AE >) ]> 

3 i 

where Nj is the yield of category j (signal, continuum background, reflection, r]K*/r]p), Pj(M^ C i, AE{) is the prob- 
ability density for the «th event and N is the total number of events. The PDFs for the signal, the reflection 
background and the r\K* jr\p feed-down are modeled with two-dimensional M^ C -AE smooth functions obtained using 
MC. The signal peak positions and resolutions in Mb c and AE are adjusted according to the data-MC differences 
using large control samples of B — > Dir and D° — > K + tt~it /tt°tt° decays. The continuum background in AE is 
described by a first or second order polynomial while the M^c distribution is parameterized by an ARGUS function, 
f(x) = xyj 1 - x 2 exp [— £(1 — x 2 )], where x is M hc divided by half of the total center of mass energy 0. Thus 
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the continuum PDF is the product of an ARGUS function and a polynomial, where £ and the coefficients of the 
polynomial are free parameters. Since B — ► r/K* branching fractions are well measured (~ 20 x 10~ 6 ) their 
feed-down to the r\K modes are fixed from MC in the likelihood fit. Since the decay B + — ► is experimentally 
poorly constrained, the amount of this background in the rjir modes is allowed to float in the fit. In the charged 
B modes, the normalizations of the reflections are fixed to expectations based on the B + — > r]K + and B + — > ?y7r + 
branching fractions and K + <-> tt + fake rates, measured using D° — > K + ir~ data. The reflection yield is first estimated 
with the assumed r/K + and 7771" + branching fractions and is then recalculated according to our measured branching 
fractions. No BB contributions are included for the B° — > r\r\ mode. 

TABLE I: Detection efficiency (e), product of daughter branching fractions ( II yield, significance (Sig.), measured 
branching fraction (B), the 90% C.L. upper limit (UL) and Acp for the B — > r\h decays. The first errors in columns 4, 6 and 
8 are statistical and the second errors are systematic. 



Mode 


e (%) n&(%) 


Yield 


Sig. 


B(10- b ) 


UL(10" b ) A C p 


B + -> rjn + 








8.0 


4.8 ±0.7 ±0.3 


0.07 ±0.15 ±0.03 




23.3 


39.4 


73.4l^;f ± 2.0 


7.0 


5.31q;°, ±0.3 


0.11 ±0.17 ±0.03 




14.8 


22.6 


19.61™ ±0-7 


3.9 


3-8111 ±0.3 


n i 1+0.35+0.04 

"■^^-o.ss-o.os 


B + -> r]K + 








3.7 


2.1 ±0.6 ±0.2 


-0.49 ±0.31 ±0.07 


r) 11 K + 


21.1 


39.4 


28.019°!° ± 16 


3.3 


2.2l°;?±0.2 


-0.45lol? ± 0.07 




13.8 


22.6 


7A±li ±0.5 


1.7 


1.5±J;J±0.2 


n 70+I.O3+O.II 
u - ,o -0.76-0.12 


B° -> V K° 








0.4 


0.3l°;?±0.1 


< 2.0 




22.9 


13.6 


-1.911;? ±°- 3 




-0.4lg;?±0.1 




mnK° 


12.2 


7.8 


3.5±l% ± 0.2 


1.3 


2-4l? j ±±0.3 




B° -» 777T 








1.8 


1.2 ±0.7 ±0.1 


< 2.5 


77 77 7T° 


17.0 


39.0 


1 o n+8.9+0.8 
lo - z -8.0-0.7 


2.4 


1.8±g | ±0.2 






11.2 


22.3 


-3.0l|;g ± 0.3 




-0.81 j;| ± 0.1 




B° -> rm 








1.1 


o.7±8:I ±o.i 


< 2.0 




16.9 


15.5 


-l.stfl ±0.1 




-0.4l°;I ±0.0 






11.3 


17.8 


7.3±t;S ± 0.2 


2.2 


2.3l?1 ±0-2 




»?3irf?37r 


7.7 


5.1 


0.31?;° ±0.1 


0.2 


0.5l?;g ±0.1 





In Table d we show the measured branching fractions for each decay mode as well as other quantities associated 
with the measurements. The efficiency for each mode is determined using MC simulation and corrected for the 
discrepancy between data and MC using the control samples. The only discrepancy we find is the performance of 
particle identification, which results in a 4.3% correction for the tjtt + mode and 1.7% for B — > i]K + . The combined 
branching fraction of the two 77 decay modes is obtained from a simultaneous likelihood fit to all the sub-samples 
with a common branching fraction. Systematic uncertainties in the fit due to the uncertainties in the signal PDFs are 
estimated by performing the fit after varying their peak positions and resolutions by one standard deviation. In the 
•qK modes, we also vary the expected r]K* feed-down by one standard deviation to check the yield difference. The 
quadratic sum of the deviations from the central value gives the systematic uncertainty in the fit, which ranges from 
3% to 6%. For each systematic check, the statistical significance is taken as the square root of the difference between 
the value of — 21n£ for zero signal yield and the best-fit value. We regard the smallest value as our significance 
incuding the systematic uncertainty. The number of B + B~ and B°B° pairs are assumed to be equal. 

The performance of the TZ requirement is studied by checking the data-MC efficiency ratio using the B + — > D°tt + 
control sample. The obtained error is 2.4-3.5%. The systematic errors on the charged track reconstruction are 
estimated to be around 1% per track using partially reconstructed D* events, and verified by comparing the ratio 
of 77 — > 7r + 7T~7r° to i] — > 77 in data with MC expectations. The 7r° and •q 11 reconstruction efficiency is verified by 
comparing the 7r° decay angular distribution with the MC prediction, and by measuring the ratio of the branching 
fractions for the two r\ decay channels: 7] — > 77 and r\ — > 7r°7r°7r . We assign 3.5% error for the ir° and ?7 77 
reconstruction. The Kg reconstruction is verified by comparing the ratio of D + — > K s tt + and D + — * K~tt + it + 
yields. The resulting Kg detection systematic error is 4.4%. The uncertainty in the number of BB events is 1%. The 
final systematic error is obtained by first summing all correlated errors linearly and then quadratically summing the 
uncorrelated errors. 

Figure^shows the M\, c and AE projections after requiring events to satisfy —0.10 GeV < AE < 0.08 GeV (—0.15 
GeV < AE < 0.10 GeV for the r7 77 and 7/7r° modes) and Afbc > 5.27 GcV/c 2 , respectively. No significant signals are 
observed for the neutral B meson modes; for these modes we set branching fraction upper limits at the 90% confidence 
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M bc (GeV/c 2 ) AE (GeV) 

FIG. 1: M bc and AE projections for (a,b) B+ -> V ty+ , (c,d) B+ -> t?A'+, (e,f) B° -> r)K° , (g,h) B° -> 777t° and (i,j) S° 7777 
decays with the ?7 77 and 773^ modes combined. Open histograms are data, solid curves are the fit functions, dashed lines show 
the continuum contributions and shaded histograms are the r\K* /r/p contributions. The small contributions around Mbc = 5.28 
GeV/c 2 and AE = ±0.05 GeV in (a)-(d) are the backgrounds from B + -> -qK + and B + r]ir + . 



level. The upper limit for each mode is determined using the combined likelihood for the two rj decay channels with 
the reconstruction efficiency reduced by la. We vary the signal PDF and the expected r\K* feed-down in the r/K 
mode to compute the likelihood as a function of branching fraction; the largest branching fraction that covers 90% of 
the likelihood area is chosen to be the upper limit. 

Significant signals arc observed for charged B decays. We investigate their partial rate asymmetries by extracting 
signal yields separately from the B + and B~ samples. Unbinncd maximum likelihood fits arc performed independently 
for the two 77 decay modes in order to reduce the systematic uncertainties. The same signal and background PDFs as 
used in the branching fraction measurement are applied. The parameters of the continuum PDF are fixed according 
to the branching fraction results. Contributions from BB backgrounds arc required to be equal for the B + and B~ 
samples. Figure |21 shows the M^ c and AE projections. The Acp results for the two 77 decay modes are combined 
assuming that the errors are Gaussian. Systematic errors due to uncertainties in the signal PDF are estimated by 
varying the peak positions and resolutions. We also check the Aqp values after varying the amount of the expected 
i]K* feed-down and the reflection background. The BB contributions are allowed to be different for the two samples 
to obtain the systematic error. The largest uncertainty is the asymmetry of the reflection. A possible detector bias in 
Acp is studied using B — ► Dtt + decays. The obtained uncertainty is 0.5%. Each Acp deviation is added quadratically 
to provide the total systematic uncertainty 

In summary, we have observed B + — > jjtt + and found evidence for B + — > rjK + ; the measured branching fractions 
and partial rate asymmetries are summarized in Table [I] We conclude that the B + — > rjTT + branching fraction is 
larger than that of B + — > rjK + . The measured B + — > rjTT + branching fraction is consistent with an earlier result 
published by the BaBar Collaboration; however, unlike the large negative Aqp measured by BaBar, the central value 



6 



in this analysis is small and positive, and is consistent with no asymmetry. For the decay B + — > rjK + , the measured 
branching fraction is 40% lower than the published result of the BaBar experiment, corresponding to a 1.3 a deviation. 
It is interesting to note that although the errors are still large, both experiments suggest a large negative Acp value 



for B + — » r\K + , which is anticipated by some theories [lif 
and upper limits at the 90% confidence level are given. 



No significant signals are found in neutral B — > r/h decays 
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FIG. 2: Afbc and AE projections for (a,b) B + -t r]TY + , (c,d) B~ -» ?77r~, (e,f) B+ -> 77^+, and (g,h) B~ -> jyif" with the 
?7 77 and 7)3^ modes combined. Open histograms are data, solid curves are the fit functions, dashed lines show the continuum 
contributions and shaded histograms are the r\K* /rjp contributions. Small curves around M bc = 5.28 GeV/c 2 and AE = ±0.05 
GeV are the reflection background on B + —* r/Tv + and B + — > r)K + . 
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